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1.  Objective 

The  objective  of  this  project  was  to  investigate  theoretically  the  microscopic  structure  and 
dynamical  behavior  of  dilute  supercritical  mixtures.  This  was  accomplished  using  two 
methods:  integral  equation  calculations  and  molecular  dynamics  simulations.  The  former 
provides  information  on  the  spatial  distribution  of  molecular  species;  the  latter  also  provides 
information  on  dynamics. 

2.  Main  Findings 

2.1  Integral  Equation  Calculations 

Two  highly  dilute  binary  supercritical  systems  were  studied:  one  attractive  and  one 
repulsive.  The  spherically-symmetric  Lennard-Jones  potential  was  used  in  both  cases. 
According  to  the  classification  introduced  by  Debenedetti  and  co-workers  (Petsche  and 
Debenedetti,  1991),  attractive  mixtures  consist  of  a  large,  non-volatile  solute  dissolved  in  a 
volatile  supercritical  solvent;  the  solute’s  characteristic  interaction  energy  is 
correspondingly  larger  than  the  solvent’s.  Repulsive  behavior  occurs  when  the  solvent  is 
larger  than  the  solute  and  has  a  correspondingly  higher  characteristic  interaction  energy.  In 
attractive  mixtures,  there  is  long-range  solvent  enrichment  around  solute  molecules  in  the 
vicinity  of  the  solvent’s  critical  point;  solubility  is  enhanced  with  respect  to  ideal-gas 
behavior  at  the  same  temperature  and  pressure;  and  the  solute’s  partial  molar  volume  and 
enthalpy  are  large  and  negative,  diverging  to  at  the  solvent’s  critical  point  in  the  limit  of 
infinite  dilution.  In  a  repulsive  mixture,  there  is  long-range  solvent  depletion  around  solute 
molecules  in  the  neighborhood  of  the  critical  point,  and  the  solute’s  partial  molar  volume 
and  enthalpy  are  large  and  positive,  diverging  to  +°°  at  the  solvent’s  critical  point  in  the 
limit  of  infinite  dilution.  Potential  parameters  for  both  mixtures  are  listed  in  tables  I  and  H 

Table  I.  Lennard-Jones  Parameters  for  Naphthalene  in  C02 

Naphthalene  in  C02  oi{  (A)  aii/q22  £8{/k  (K)  e,i/£22 

solute-solute  j  11  6.199  1.634  554.4  2.458 

solvent-solvent  22  3.794  1.000  225.5  1.000 

solute-solvent  RZ^ . 4.996  1..1I7 . ~  3534“ .  4.5677 

Table  II.  Lennard-Jones  Parameters  for  Neon  in  Xenon 


kjviuw  uviwvv  ,  i  w  v  V.  :7T.; . 

solvent-solvent  22  4.047  1.000  231.0  1.000 

solute-solvent  f  12  3.433  0.848  87.0  0.377 


Neon  in  Xenon  ij 

cr»1nfy=>_cnlntt»  1  1 


U 


9  890 


°ij/a22 


ei 


0  697 


£u(£22 


39;8 


0.142 
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The  mixtures  listed  in  Tables  I  and  II  have  been  called  naphthalene  (solute)  in  carbon 
dioxide  (solvent)  (Cochran  and  Lee,  1989)  and  neon  (solute)  in  xenon  (solvent)  (Petsche 
and  Debenedetti,  1989).  We  use  this  terminology,  although  it  should  be  understood  that 
except  for  yielding  reasonable  estimates  for  pure-component  critical  temperatures  and 
densities  with  the  parameters  listed  in  Tables  I  and  II,  the  Lennard-Jones  potential  is  a 
rudimentary  representation  of  the  actual  interaction  potentials. 

The  most  important  findings  that  resulted  from  this  work  were: 

(a)  Large  enhancements  in  the  calculated  local  density  of  solute  molecules 
around  a  central  solute  molecule.  Local  densities  roughly  an  order  of  magnitude  in 
excess  of  the  bulk  density  were  found.  This  is  illustrated  in  Figure  1. 

(b)  Local-bulk  asymmetry  is  most  pronounced  at  sub-critical  densities.  The 

maximum  enhancement  in  the  solute-solute  distribution  function  shown  in  Figure  1  occurs 
at  a  reduced  density  of  0.58.  The  fact  that  such  pronounced  local-bulk  asymmetries  are 
more  important  at  sub-critical  densities  had  not  been  previously  understood. 

(c)  The  solute’s  chemical  potential  is  determined  by  short-ranged  solvation 
interactions.  This  is  illustrated  in  Figure  2,  which  shows  the  cutoff  distance  (beyond 
which  the  mixture  is  assumed  to  be  homogeneous)  needed  to  obtain  specified  percentages 
of  the  asymptotic  value  of  the  solute’s  fugacity  coefficient. 

(d)  Large  depletions  in  the  distribution  of  solvent  molecules  around  the 
solute  exist  in  repulsive  mixtures.  This  is  illustrated  in  Figure  3. 

Taken  together,  (a)  and  (b)  imply  that  solvation  and  local-bulk  asymmetry  are  unrelated  to 
proximity  to  critical  points.  These  short-ranged  effects  are  a  consequence  of  differences  in 
size  and  characteristic  interaction  energy  between  the  solute  and  the  solvent. 


2.2  Molecular  Dynamics  Simulations 

The  objective  of  this  study  was  to  understand  the  mechanisms  of  solvation  in  supercritical 
solvents,  and  in  particular  to  understand  in  what  ways  the  translational  and  rotational 
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motion  of  solutes  is  affected  by  high  local  solvent  densities.  The  system  benzene  (solute)  - 
carbon  dioxide  (supercritical  solvent)  was  investigated.  Both  molecules  were  modeled  as 
rigid  polyatomics  interacting  via  multi-site  Lennard-Jones  potentials.  In  addition, 
interactions  between  carbon  dioxide  molecules  included  a  point  quadrupole  moment.  Site 
parameters  are  listed  in  Table  III. 


Table  III.  Site 

(a) 

Simulations 

Potential 

Parameters  for 

Benzene-Carbon  Dioxide 

Molecule 

£/k  (K) 

a  (A) 

Carbon  Dioxide  f 

C-C 

29 

2.785 

:  ■,  :  J 

Benzene 

0-0 

83.1 

:  .  .iyyy... 

3.014 

:  f 

CH-CH 

95.05 

3.350 

(a)  Lorentz-Berthelot  used  for  unlike  pairs.  CO  distance  =  1.16  A;  CH-CH  distance  — 
1.765  A 


Simulations  were  done  at  two  reduced  temperatures  (T/Tc  =  1.05,  1.1)  and  five  reduced 
densities  (p/pc  =  0.5,  0.75,  1,  1.5,  2).  In  both  cases,  the  critical  parameters  are  the 
solvent’s. 

The  main  results  of  this  work  were: 

(a)  The  distribution  of  solvent  molecules  around  benzene  is  highly 
asymmetric.  Benzene  is  preferentially  solvated  in  the  polar  position,  where  the  solvent 
density  can  be  as  high  as  seven  times  its  bulk  value.  This  is  shown  in  Figure  4.  Note  that 
the  difference  between  local  and  bulk  conditions  is  more  pronounced  at  the  lowest  density 
investigated.  No  unusual  solvation  effects  are  seen  at  the  critical  density. 

(b)  Benzene’s  translational  diffusion  in  supercritical  carbon  dioxide  scales 
inversely  with  bulk  density,  and  is  unaffected  by  proximity  to  the  critical 
point.  This  is  shown  in  Figure  5. 

(c)  The  rotational  motion  of  benzene  is  remarkably  insensitive  to  bulk 
density  over  the  range  of  densities  explored.  No  transport  anomalies  were 
found  near  the  solvent’s  critical  point. 
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(d)  The  diffusive  exchange  of  solvation  shell  molecules  is  described  by  a 
single  exponential,  with  a  time  constant  that  increases  monotonically  with 
increasing  density.  This  time  constant  ranges  from  3  psec  at  50%  of  the 
critical  density  to  7  psec  at  twice  the  critical  density.  This  is  shown  in  Figure  6. 
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Figure  1:  Integral  equation  calculation  of  the  density  dependence  of  the  solute-solute  pair 

correlation  function  for  the  naphthalene-carbon  dioxide  system.  The  vertical  axis  is  the  ratio 

of  local  to  bulk  concentration  of  naphthalene  molecules  around  a  central  naphthalene 

molecule.  The  transverse  axis  is  the  distance  from  the  central  naphthalene  molecule,  in  units 

* 

of  the  Lennard-Jones  diameter  of  carbon  dioxide  (Table  I).  p  is  the  number  density 

(molecules  per  unit  volume),  in  units  of  the  cube  of  the  Lennard-Jones  diameter  of  carbon 

* 

dioxide.  The  critical  density  of  the  pure  solvent  is  at  p  =0.31.  The  temperature  is  7% 

-9 

above  the  solvent’s  critical  point,  and  the  solute  mole  fraction  is  10  .  Note  the  large  first 
peak,  especially  at  sub-critical  densities.  (Tom  and  Debenedetti,  1993). 
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Figure  2:  Integral  equation  calculation  of  the  cutoff  distance  needed  to  obtain  the  specified 
percentages  of  the  asymptotic  value  of  the  solute’s  fugacity  coefficient.  System  naphthalene 
(solute)-carbon  dioxide  (solvent).  The  cutoff  (vertical  axis)  is  in  units  of  carbon  dioxide’s 
Lennard-Jones  diameter  (Table  I).  Temperature  and  mole  fraction  as  in  Figure  1.  Density 
units  as  in  Figure  1.  The  cutoff  distance  calculated  here  is  the  solvation  region  that  defines 
the  solute’s  chemical  potential.  Thus,  the  calculations  show  that  solubility  is  determined  by 
the  local  environment  within  1 1  to  20  A  around  a  central  naphthalene  molecule.  (Tom  and 
Debenedetti,  1993). 
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Figure  3:  Integral  equation  calculation  of  the  density  dependence  of  the  ratio  of  local-to- 
bulk  density  of  supercritical  xenon  atoms  around  a  central  neon  atom.  The  density  is 
averaged  over  a  sphere  of  specified  radius,  expressed  in  units  of  the  Lennard-Jones 

diameter  of  xenon  (Table  II).  p  is  the  number  density  (molecules  per  unit  volume),  in 

units  of  the  cube  of  the  Lennard-Jones  diameter  of  xenon.  The  critical  density  of  the  pure 
* 

solvent  is  at  p  =  0.31.  The  temperature  is  7%  above  the  solvent’s  critical  point,  and  the 

.9 

solute  mole  fraction  is  10  .  Note  the  large  depletion  region  where  the  local  density  is 
smaller  than  the  bulk.  This  is  due  to  the  mixture’s  repulsive  behavior.  (Tom  and 
Debenedetti,  1993). 
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Figure  4:  Sector-dependent  radial  distribution  function  of  supercritical  carbon  dioxide  around  benzene,  g  is 
the  ratio  of  local  to  bulk  density.  Also  shown  is  the  geometry  for  the  definition  of  sectors.  T*  =  kT/ecc;  p‘  = 
pa3cc  (see  Table  m  for  ecc  and  acc  values).  T  (K)  =  29T*;  p  (mol/liter)  =  76.87  p*.  For  the  potential 
parameters  used  in  this  study,  Tc  =  330.6  K  (T*  =  11.4);  pc  =  11.15  mol/liter  (p  =  0.145).  The  two  salient 
features  are  the  preferential  solvation  by  carbon  dioxide  in  the  polar  region,  and  the  fact  that  the  most 
pronounced  enhancements  in  local  solvent  density  occur  at  sub-critical  densities.  Molecular  dynamics 
simulation.  (Inomata  et  al.,  1996). 


Figure  5:  Translational  diffusion  coefficient  of  benzene  in  carbon  dioxide  at  two 
supercritical  temperatures.  Vertical  arrow  indicates  carbon  dioxide’s  critical  density. 

Experimental  data  are  from  Bueno  et  al.,  J .Chem. Eng. Data,  38,  344  (1993).  T  =  kT/ecc; 

p  =  pa3cc  (see  Table  III  for  ecc  and  Gcc  values).  T  (K)  =  29T  ;  p  (mol/liter)  =  76.87  p  . 

* 

For  the  potential  parameters  used  in  this  study,  Tc  =  330.6  K  (T  =  11.4);  pc  =  11.15 
* 

mol/liter  (p  =  0.145).  Note  the  smooth  variation  of  the  translational  diffusivity  in  the 
critical  region.  Molecular  dynamics  simulation.  (Inomata  et  al.,  1996). 
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Figure  6:  Decay  of  the  fraction  of  tagged  carbon  dioxide  molecules  within  a  sphere*  of 
radius  6.7  A  centered  at  benzene’s  center  of  mass  that  remain  in  the  sphere  at  time  t.  T  = 

kT/ecc;  p*  =  pa3cc  (see  Table  III  for  ecc  and  cxcc  values).  T  (K)  =  29T  ;  p  (mol/liter)  = 
*  * 

76.87  p  .  For  the  potential  parameters  used  in  this  study,  Tc  =  330.6  K  (T  =  1 1.4);  pc  = 
* 

11.15  mol/liter  (p  =  0.145).  The  decay  is  described  by  a  single  exponential.  Note  the 
slowing  down  of  the  decay  rate  upon  increasing  the  density. 
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